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Abstract 

The pion and kaon electromagnetic form factors Fm{Q'^) are cal- 
culated at the leading order of pQCD using the running coupling con- 
stant method. In computations dependence of the mesons distribu- 
tion amplitudes on the hard scale is taken into account. The Borel 
transform and resummed expression for Fm{Q'^) are found. The effect 
of the next-to-leading order term in expansion of as{XQ'^) in terms 
of Q!s(Q^) on the pion form factor Ft^{Q'^) is discussed, comparison is 
made with the infrared matching scheme's result. 



*Talk given at the Euroconference QCD 98, Montpellier 2-8th July 1998, France, to 
appear in Proceedings. 
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1. Investigation of the infrared (ir) renormalon effects in various inclusive 
and exclusive processes is one of the important and interesting problems in 
the pQCD [|I| (and references therein). It is known that all-order resumma- 
tion of ir renormalons corresponds to the calculation of the one-loop Feyn- 
man diagrams with the running coupling constant as{—k'^) at the vertices 
or, alternatively, to calculation of the same diagrams with non-zero gluon 
mass. Both these approaches are generalization of the Brodsky, Lepage and 
Mackenzie (BLM) scale-setting method and amount to absorbing certain 
vacuum polarization corrections appearing at higher-order calculations into 
the one-loop QCD coupling constant. Studies of ir renormalon problems 
have opened also new prospects for evaluation of higher twist corrections to 
processes' characteristics. 

Unlike inclusive processes exclusive ones have additional source of ir renor- 
malon contributions 0-^; integration over longitudinal fractional momenta 
of hadron constituents in the expression of the elm form factor generates such 
corrections. 

2. In the context of pQCD the meson elm form factor has the form, 
Fm{Q') = C t dxdycj)l,{y, Q')Th{x, y; Q\ as(Q'))0M(a;, Q'), (1) 
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where = —q^ is the square of the virtual photon's four-momentum. Here 
(f)M is the meson distribution amplitude, containing all non-perturbative 
hadronic binding effects, whereas Th is the hard-scattering amplitude of the 
subprocess qq' + 7* ^ qq' and can be found using pQCD. In (1) Q"^ is the 
factorization and renormalization scale, which is taken as the square of the 
momentum transfer of the exchanged hard gluon in corresponding Feynman 
diagrams. Such choice for allows one to remove large terms proportional 
to ln(Q^//i|.) and ln((5^/yu|.) from the next-to-leading order correction to Th 

i- 

At the leading order Th is given by the following expression 



Th- 
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One of the important ingredients of our study is the choice of the meson 
distribution amplitude (pAiiXjQ"^). In the framework of pQCD it is possible 
to predict the dependence of 4>m{x, Q^) on using evolution equation, but 
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not its shape (its dependence on x). In this work for the pion and kaon we 
use model distribution amphtudes proposed in Refs.0-0. For the pion they 
have the following form 



M^, nl) = (f)l,y{x) [a + h{2x - if + c(2x - 1)' 
where (plsyi.^) is the pion asymptotic distribution amplitude, 



(3) 



(4) 



and /tt = 0.093 GeV is the pion decay constant. 

The constants a, 6, c in (3) were found by means of the QCD sum rules 
method at the normalization point /iq = 0.5 GeV and take different values 
0] 5 HI 5 defining the pion's alternative distribution amplitudes. For the kaon 
we have 



/io) = + - 1) + c(2x - 1)' 



(5) 



with a = 0.4, 6 = 3, c = 1.25 and /k = 0.122 GeV. 

The dependence of the meson distribution amplitude on can be ob- 
tained by means of the following expression 



'x)Y^rr,Cl/\2x-l)Ar^ 



(6) 



n=0 



Here |C^/^(2a; — 1)} are the Gegenbauer polynomials, 7„ is the anomalous 
dimension and 



as{Q': 
as(/io) 



(3q = 11 — 2nj/3 being the QCD beta-function's one-loop coefficient. 

The QCD running coupling constant cxsiQ'^) in Eq-(2) suffers from ir 
singularities associated with the behaviour of the otsiQ"^) in the soft regions 
X— ^0, y^O; X ^ 1, y ^ 1. Therefore, Fm{Q'^) can be found after proper 
regularization of otsiQ"^) in these soft end-point regions. In the framework of 
the frozen coupling approximation such regularization is achieved by equating 

to its mean value (5^/4 and removing asiQ"^/^) as the constant factor 
from integrand in (1). As a result, one obtains form factors with the same 



3 



shape, but different magnitudes depending on a distribution amphtude used 
in computations. To solve this problem in the context of the running coupling 
method let us relate the running coupling asi^XQ"^) in terms of OisiQ"^) by 
means of the renormalization group equation [[10] 



as 



al l3i ln[l + (a5/3o/47r)lnA] 

(«5/?o/47r)lnA An /?o[l + (a5/?o/47r) In A]^ ' 



(7) 
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where as is the one-loop QCD coupling constant asiQ"^) and /?i 
38?2j/3 is the beta-function's two- loop coefficient. 

3. As was demonstrated in our works integration in (1) using 

(2) and (7) generates ir divergences and as a result for Fm(Q^) we get a 
perturbative series with factorially growing coefficients. This series can be 



resummed using the Borel transformation |TT 



(16vr/M) 



/5o 



rfwexp 
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(8) 



where B[Q'^Fm]{u) is the Borel transform of the corresponding perturbative 
series §, §. 

When we take both of the variables x,y in Eq.(2) as the running ones, 
for B[Q^FM]iu) we find 



B 



Q^Fm 



N 

E 

n=l 



(n — u) 



+ 



n — u 



(9) 



The exact expressions for m„, 1„ can be found in Ref . |T2[ . 

The Borel transform (9) has double and single poles aX u = n. They are 
ir renormalon poles, which are responsible for divergence of the perturbative 
series for Q'^Fm{Q'^)- The resummed expression (8) can be calculated with 
the help of the principal value prescription |lTO|,[p!T[| 



Q^Fm{Q' 



;i67r/ 



M) 



\2 N 



E 

n=l 



h (1„ + m„lnA)— — 

n A" 



where li{X) is the logarithmic integral ||T3 

h{X) = P.V. f'-^, X = Q^/K\ 
Jo mx 



(10) 



(11) 
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In (10) we have taken into account the dependence of the distribution 
amphtude 4>m{x, Q"^) on the scale Q"^ (in [12], for simphcity, we have replaced 
(Pm{x,Q'^) (Pm{x,Q'^)) and it is the generalization of our results for the 
pion and kaon elm form factors [4] , [5] . 

In the case with one frozen (for example, y) and one running (x) variables, 
which corresponds to the choice = Q^{1 — x)/2 and = Q'^x/2 in (2), 
for the pion we find (asymptotic distribution amplitude (4)) 



2Po 



X A2 



A = gV2A2. (12) 



The ir renormalon analysis carried out (Eqs.(10),(12)) allows one to estimate 
power corrections to the light mesons' form factor. Another way of such 
estimation is the ir matching scheme, in the context of which, one explicitly 
divides power corrections from a full expression by introducing moments of 
as at low scales as new non-perturbative parameters fl^ . As an example, let 



us consider the pion form factor {4>asyi^))- By freezing one of variables (y) 
we can express Q^F^^ in terms of moment integrals F„ defined by the formula 



Qp Jo 

After simple calculations we get 



4 



+^ [1 - 2r(l, 2z) + r(l, Az)] + ^ [3 - 4r(2, 2z) + r(2, 4^)]| , (13) 
where fi = 2 GeV is the ir matching scale, z = \ia{Q/ ^/2|J?j, asiQ'^/2) and 



r(n, x) is the incomplete gamma function [|T^]. 

4. Results of our numerical calculations are shown in Fig.l. For com- 
putation of Q^Ftt in the context of the ir matching scheme, we use N = 4 
perturbative estimate for moment integrals Fp[Q) (in Eq.(13), expression for 
= 1 is written down). 

It is worth noting that non-perturbative parameters F2{fi) and F^lfi) can 
be found using the running coupling method 



5 



Po A f3o 

For A = 0.2 GeV we get F2(/i = 2GeV) = 0.421, F^ifx = 2GeV) = 0.347. 
In our calculations we take for F2{fi = 2GeV) ~ 0.5 the value deduced 
from experimental data and for F4-running coupling prediction 0.347. As is 
seen, Eq.(12) and ir matching scheme give approximately the same results, 
excluding a region of small Q^. 

In Fig. 2, as an example, the dependence of the ratio Rk = [Q'^Fk{Q'^)Y'^'^ / 
on Q^, where [(5^-FV(<5^)]° is the kaon form factor in the frozen 
coupling approximation, is shown. Here we take into account a dependence 
of the kaon's distribution amplitude on (curve 1, the dashed curve is 
taken from Ref.[5], where this dependence was neglected). The ir renormalon 
corrections can be transferred into the scale of asiQ^) in [Q^Fk{Q^)]^ 

Q2 ^ ^m')Q2^ f^Q2^ = + C2as{Q') + C3«|(g'). (M) 

Numerical fitting allows us to get (forrij = 3) Ci ~ —1.304, C2 — —35.604, C3 ~ 
127.25 (curve 2). The similar results can be obtained also for the pion. 

It is interesting to clarify an importance of the second term in (7) for 
considering problem. Investigations demonstrate that an effect of the second 



term on a whole result is small (see, Ref. |jT5| for details). 

5. It is evident that ir renormalon effects enhance the ordinary (frozen 
coupling) perturbative predictions for the pion, kaon elm form factors approx- 
imately two times. Our recent studies confirm that neither a dependence of 
4^m{x,Q'^) on nor the next-to-leading order term in Eq.(7) changes this 
picture considerably. This feature of ir renormalons may help one in solution 
of a contradiction between theoretical interpretations of experimental data 
for the photon-to-pion transition form factor F^^r from one side and for the 
pion elm form factor F^^ from another side. Thus in Ref. |l^ the authors noted 



that the scaling and normalization of F^^^ to favor of the pion asymptotic-like 
distribution amplitude. But then prediction for obtained using the same 
amplitude is lower than the data by approximately a factor of 2. We think 
that in this discussion a crucial point is a chosen method of integration in 
(1). Indeed, unlike F.^ the expression for F^tt at the leading order of pQCD 
does not contain an integration over as{xQ'^). In other words, the running 
coupling method being applied to Eq.(l) and to the expression for (see. 
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Ref. ||16||) enhances the perturbative result for the pion elm form factor and 
at the same time, does not change F^t,. This allows us to suppose that 
in the context of pQCD the same pion distribution amplitude may explain 
experimental data for both F^^^ and F^^. 
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FIGURE CAPTIONS 

Fig.l The pion elm form factor calculated using (f)lsy{x) (4). The curves 
correspond to the following computational schemes: 1- resummed expression 
(10), 2-resummed expression (12), dashed-ir matching scheme, dot-dashed- 
frozen coupling approximation. 

Fig. 2 The ratio Rk as a function of Q^. 
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